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^ , Abstract: 

^ ' Continuous phase transitions in spin systems can be formulated as percolation of 

suitably defined clusters. We review this equivalence and then discuss how in a similar way, 
the color deconfinement transition in SU(2) gauge theory can be treated as a percolation 
phenomenon. In the presence of an external field, spin systems cease to show thermal 
critical behavior, but the geometric percolation transition persists (Kertesz line). For 
if 7^ 0, we study the relation between percolation and pseudocritical behavior, both for 
continuous and first order transitions, and show that it leads to the necessity of an in- 
dependent cluster definition. A viable formulation of this kind could serve as definition 
^ ! of deconfinement in QCD with dynamical quarks. 
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1. Introduction 



Critical phenomena in spin systems and elsewhere are usually described in terms of 
spontaneous symmetry breaking or of singular behavior of the partition function. On the 
other hand, the size of domains of parallel spins generally increases near criticality, and so 
the search for a geometric "connectivity" formulation has been a long standing challenge. 
For some classes of systems, the problem is solved, and I begin by recalling how this is 
. done. 

rS I The Ising model is defined by the Hamiltonian 



niT,H) = -JY.s,s,-HJ2s^■, Si = ±l (1) 

<ij> i 

which is used to calculate the partition function 

Z(T,ii) = ^exp{-/57^(T,ii)}; P ^ 1/T; (2) 

here J denotes the nearest neighbor spin-spin coupling and H an external field. Given 
Z[T,H), the magnetization and (isothermal) susceptibility are obtained as 

m{T,H) = {d log Z/dH)T, x{T,H) = {dm/dH)T. (3) 
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For vanishing external field, H = 0, the Ising model leads to a critical point T = T^. 
for T < Tc, the Z2 symmetry of the Hamiltonian 7i(T, 0) is spontaneously broken by the 
state of the system. The magnetization becomes the order parameter for the transition 
and vanishes for T Tc, while the susceptibility diverges there: 

m(T,0)~(r,-Tf, T<T,; 0) ~ |T, - T^. (4) 

The critical exponents f3 and 7 specify the critical behavior and define the universality 
class of the transition. 

It was shown by Fortuin and Kasteleyn that the partition function of the model, 
given in Eq. (0) as sum over spin configurations, can for H = he equivalently written 
as a sum over configurations of clusters containing n lattice points. In the resulting form 

nij=l nij=0 

Z{T,0)=J:{ Up^U (1-P.)2^(")}, (5) 

n <ij> <ij> 

Pb = 1 — exp{— 2/3J} defines a bond weight between adjacent sites, while riij = 1(0) 
specifies whether two sites are bonded or not. The number of ra-clusters is given by c{n). 

In percolation theory, the observable corresponding to the order parameter is the 
percolation strength P(T), defined as the probability that a given site belongs to a per- 
colating cluster. In place of the susceptibility we now have the average cluster size S{T); 
in the percolation region only finite clusters are included in the averaging. Using the FK 
formalism, Coniglio and Klein showed that the percolation point coincides with the 
thermal critical temperature T^, and that the percolation behavior is governed there by 

P(T) ~ (T, - Tf, T < T,; 5(T) ~ |T, - T\\ (6) 

where (3 and 7 are the usual thermal Ising exponents. For the Ising model, thermal critical 
behavior can thus be equivalently described as percolation of FK clusters. 

A non-vanishing external field, H ^ 0, explicitly breaks the previous Z2 symmetry 
of the Hamiltonian, and as a result, partition function and thermodynamic observables 
become analytic: there is no more thermal critical behavior. The magnetization now 
never vanishes; H aligns some spins at any temperature (see Fig. |I|). 

However, the percolation transition persists also for H 0. For any value of H 
there exists a percolation point Tp{H), such that for T < Tp there is percolation and for 
T > Tp not 0. The critical behavior along the "Kertesz line" Tp{H) (see Fig. |l|) is now 
determined by percolation critical exponents, so that the transition is no longer in the 
Ising universality class. 

We note here that the equivalence between thermal and percolation behavior can in 
fact be restored if clusters are not required to be spatially connected. The introduction of 
an if-dependent ghost spin not on the lattice but able to couple to any lattice site allows 
spatially disjoint sites to form a cluster ||^ and thus leads to percolation for all if 7^ 0. 
We are here interested in the relation between spatial connectivity and thermal behavior; 
therefore we shall not pursue the ghost spin approach any further. 

One major reason for our study is to address critical behavior in finite temperature 
QCD; let us therefore elaborate a little the problems which arise there. The constituents 
of strongly interacting matter, the hadrons, are finite size color-neutral bound states of 
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Figure 1: Magnetization and Kertesz line in the Ising model. 

two (mesons) or three (nucleons) colored quarks coupled by colored gluons. At high den- 
sity, when the hadrons strongly overlap, each quark finds in its immediate vicinity many 
others, so that a partitioning into subsets of two or three quarks in a given volume no 
longer makes any sense, suggesting that we now have quark matter. High temperatures 
lead through particle-antiparticlc production to high density. Hence we expect at suffi- 
ciently high temperatures a transition from hadronic matter, with color-neutral hadrons 
as constituents, to a plasma of unbound colored quarks and gluons. One of the basic 
issues of strong interaction thermodynamics is this quark-hadron transition. 

The theory of strong interactions, quantum chromodynamics (QCD), is defined by the 
QCD Lagrangian jC(mq), which contains the bare quark mass as an open parameter. In 
QCD thermodynamics, two limits are of particular interest. For massless quarks, rrig = 0, 
>C(0) is invariant under chiral transformations. This chiral symmetry is spontaneously 
broken for T < T^: through a surrounding gluon cloud, the quarks acquire a dynami- 
cally generated 'constituent quark' mass Mg, which essentially determines the mass of the 
conventional hadrons. Thus the mass of a nuclcon as three quark bound state is approx- 
imately 3Mq. At T = T^, chiral symmetry is restored, Mg and hadrons dissolve into 
quarks and gluons. In the opposite limit, for nig oo, the quarks drop out and QCD 
becomes a (non-Abelian) SU(3) gauge theory, which describes a system of interacting 
gluons. At low temperatures, these bind to form 'glueballs' as the hadrons of the theory. 
The resulting Lagrangian has a global symmetry. For some T > T^, this symmetry is 
spontaneously broken; color screening causes the colorless glueballs to dissolve into col- 
ored gluons, corresponding to the onset of color deconfinement. QCD thus leads to two 
genuine phase transitions: 

• for niq — 0, chiral symmetry restoration at T^, corresponding to the vanishing of 
the constituent quark mass; 

• for niq — oo, spontaneous breaking at Tc corresponding to the onset of color 
deconfinement in pure gauge theory. 

In the real world, ruq is small but finite; what happens then? For three massless quark 
fiavors, the Z^ transition is first order, so a discontinuity will persist for some range of 
masses m^"*^ < mq. For rUq < m®"*^, there is no further thermal critical behavior, which 
means that there is no true phase transition separating hadronic matter and quark-gluon 
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plasma (see Fig. Is the transition from one state to the other then just a rapid cross- 
over? We would like to speculate that it is the Kertesz line of QCD, corresponding to the 
onset of percolation for clusters of deconfined matter 0, and we have carried out some 
preliminary studies to address this idea. 
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Figure 2: The phase structure of QCD (left) and of the 3-d 3-state Potts model (right) 

2. Percolation and Deconfinement 

As a simplified model, we consider here pure SU(2) gauge theory, which describes a 
system of interacting gluons of three different colors. Given such a medium at temper- 
ature T, we denote by V{r,T) the potential between two static color charges Q and Q 
separated by a distance r. Below a certain critical temperature Tc, there is confinement: 
the potential rises linearly with distance, so that V{r,T) oo for r —>■ oo. Above Tc, 
color screening limits the range of the potential, so that now V{r,T) —>■ g{T) for r oo, 
with a finite g{T) > 0. The Polyakov loop expectation value, defined by 

L(T) = ;imexp{-V(r,T)/T}, (7) 

thus vanishes for T < Tc, i.e., in the confinement region, and takes on a non-zero value 
for T > Tc, when deconfinement set in. Hence it is the deconfinement order parameter. 

In the lattice formulation of finite temperature QCD, one studies systems in the usual 
three space dimensions, with one additional dimension specifying the temperature. A 
Polyakov loop becomes in lattice SU(2) an up or down spin of arbitrary amplitude, as- 
sociated to a given spatial lattice site and oriented in the temperature direction; the 
expectation value L{T) is the average over all spatial lattice sites. Polyakov loop configu- 
rations on a given spatial lattice thus form something like a continuous spin Ising model, 
and it was therefore conjectured that the deconfinement transition in SU(2) gauge theory 
is in the universality class of the Ising model 0. This conjecture was subsequently con- 
firmed by extensive lattice studies, showing that the critical exponents of the two models 
agree with remarkable precision ||^. 

We now would like to check if the deconfinement transition can also be described as 
percolation of like-sign Polyakov loop clusters. The first question to be addressed is the 
generalization of the global FK bond weight 1 — exp{— 2J/T} of the conventional Ising 
model with Sj = ±1 to the case where the amplitude of the spin becomes a local variable. 
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It was shown in both analytic and numerical studies of the continuous spin Ising model 



rO| that percolation and the thermal transition in fact become equivalent if the bond 
weights are taken to have the local form 1 — exp{ — (2J/T)sjSj} for adjacent like-sign 
spins. 

Turning to SU(2) gauge theory, we therefore start with the local FK bond weight 
p{i, j) = l—exp{2KLiLj } expressed in terms of the local values Li, Lj of like-sign Polyakov 
loops on adjacent spatial lattice sites the variable k must take into account the tem- 
perature dependent Polyakov loop coupling. The generalization to SU(2) thus leads to 
two basic problems: how can one determine the functional form of k{T), and is it justified 
to approximate the correlations provided by the theory to only nearest neighbor inter- 
actions? In the strong coupling limit [|ll[], one finds for Nj- = 2 that k{T) ~ (/3(T)/2)^, 



where (3{T) is the known temperature dependent coupling in the SU(2) gauge theory ac- 
tion. Using this relation, we thus generate SU(2) configurations at different temperatures 
and study the clustering properties. The results ^ show overall agreement, but some 
discrepancies in details; in particular, the critical couplings for the thermal transition and 
for deconfinement are obtained with very high precision and disagree slightly. 

It therefore seems worthwhile to consider a formulation in which the strong coupling 
approximation and the restriction to nearest neighbor interactions are removed. This can 
be done in an effective spin description, in which the thermal system is parametrized in 
a spin model allowing longer range interactions as well, with up to 20 different couplings 
between different spins [Q. The values of these couplings are then determined numerically 
near the thermal critical point which in turn provides generalized bond weights k„, 

p{i,i + 1) = 1 - exp{-2KiLiLi+i}, p{i,i + 2) = 1 - exp{-2KiLiLi+2}, ■■■ (8) 

The resulting agreement between the critical thermal and percolation parameters is ex- 
cellent, as seen in Table 1. The approach used here does not solve the question of how 
to define in general the FK bond weights in terms of the SU(2) Lagrangian; but it does 
show that the equivalence between thermal and geometric descriptions can be established 
numerically. 





Critical point 
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7/z/ 


V 


Percolation 


1.8747(2) 


0.528(15) 


1.985(13) 


0.632(11) 


Spontaneous Z2 Breaking 


1.8735(4) 


0.523(12) 


1.953(18) 


0.630(14) 


Ising Model |12 




0.518(7) 


1.970(11) 


0.6289(8) 



Table 1: Parameters for percolation and spontaneous Z2 breaking in 3+1 SU(2). 

We therefore conclude that color deconfinement in SU(2) gauge theory can be con- 
sidered either as spontaneous symmetry breaking or as the onset of Polyakov loop 
percolation. Both phenomena occur at the same critical temperature T^. and are described 
in terms of the same set of critical exponents. 

3. Percolation in Non-Critical Thermal Systems 

We would like to formulate a percolation description of deconfinement in full QCD with 
dynamical quarks; again we start with a preliminary study and consider spin systems with 
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a non- vanishing external field H ^ 0. The basic issue here is to generalize the FK bond 
weights to include a dependence on H: Pb(T, H). This problem is not yet solved; we shall 
address two closely related questions. 

Consider the 2-d Ising model with H 0; the behavior of the magnetization was 
shown in Fig. |l[ Included there is the Kertesz line Tp[H), separating the percolating 
(T < Tp{H)) and the non-percolating region (T > Tp{H)). In addition, one can also 
consider the pseudocritical line T^{H), defined by the temperature value for which the 
isothermal susceptibility xt = {dm/ dH)^ peaks for a given if; for if = 0, it would 
diverge at Tc. 

We want to compare these two temperature lines |]T3; in all percolation calculations, 



we use the conventional (ii-independent) FK bond weights introduced in section 1. In 
terms of the reduced variables t = (T — Tc)/Tc and h = H/ J we find for small h 

t^ = C^h^/PS tp = C2/l^ (9) 

for the pseudocritical Ising line and the Kertesz line, respectively, with constants Ci and 
C2- For the Ising model, one has the exact result 1/(36 = 8/15 ~ 0.533, while a numerical 
study of percolation gives x = 0.534(3). For small h, the two lines thus have the same 
/i-dependence. However, ci ^ C2, so that the lines do not coincide. For large h, this 
is obvious: tp{h) —> % for h oo, where U is the bond percolation temperature, but 

—>■ oo foT h oo, although the susceptibility peak becomes arbitrarily weak in this 
limit. This illustrates the necessity to obtain field-dependent bond weights Pb{T, H) for a 
reasonable description at ii 7^ 0. We note here that a comparison with the pseudocritical 
behavior of the isothermal susceptibihty is not a unique way to specify p(T,H), since 
other susceptibilities or the correlation length lead to slightly shifted peak temperatures. 

In a second study [l^, we consider the 3-d 3-state Potts model, which shows a first 



order phase transition for H = 0. When an external field is turned on, the discontinuity 
in e.g. the magnetization decreases but continues for < ii < iio, defining a line Tdis(ii). 
This ends for some value H = Hq aX Tc{Hq) with a second order phase transition which 
appears to be in the universality class of the Ising model |T^. For H > Ho, the partition 



function is analytic and there is no more thermal critical behavior (see Fig. Using the 
conventional FK bond weights p{T) = 1 — exp{— 2J/T}, we now study the percolation 
behavior of the model. The result is 

• for < ii < iio, there is a discontinuous onset of percolation at the same Tdis(ii) 
at which there is discontinuous thermal behavior; 

• at H = Hf), the percolation and the thermal transition become continuous and occur 
at the same temperature Tc{Hq). In contrast to the Ising-like thermal behavior, 
the percolation transition leads to exponents which are neither in the Ising nor in 
the random percolation universality class. This is presumably due to the infinite 
correlations at Tc{Hq), similar to what is found for pure site percolation in the 2-d 



Ising model at the Curie point |16 



Thus again there seems to exist a definite relation between percolation and thermal critical 
behavior, but with the conventional FK weight, there is not real equivalence; a bond 
weight p{T, H) seems to be needed. Here, however, in contrast to the pseudocritical 
line(s), the problem is well-defined: what p{T,H) will put percolation at Tc{Hq) into 
the Ising universality class? This is quite similar to the original question of what p{T) 
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would change random percolation enough to map the new clustering pattern onto thermal 
critical behavior, solved by the FK bond weights. If a correct p{T, H) can be determined, 
one would have full equivalence between percolation and thermal critical behavior both 
in the first order region and for the continuous transition at the endpoint; beyond that, 
percolation only would persist. Given such a framework, one could then try to address 
the similar structure in full QCD (Fig. |^). 

4. Outlook 

From what we have seen, thermal critical behavior can in many cases be formulated as 
cluster percolation. The interesting feature is, however, that percolation can persist even 
when thermal critical behavior stops, such as in spin systems with external field. Here the 
percolation strength P{T) and the cluster size S{T) remain singular when the partition 
function Z{T) becomes analytic. This raises some intriguing questions. Can one give 
a more general definition of critical behavior, such that percolation is included? What 
physical features distinguish the percolating from the non-percolating phase |0? Is it 
possible to define a more 'detailed' partition function Z{T,x) at some fixed connectivity 
measure x, which results in the usual Z{T) after summation or integration over all x? It 
seems that all these points could eventually provide a new opening to thermodynamics. 
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